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ABSTRACT 

Physical conditions in the atmospheres of tidally locked, slowly rotating hot 
Jupiters correspond to dynamical circulation regimes with Rhines scales and 
Rossby deformation radii comparable to the planetary radii. Consequently, the 
large spatial scales of moving atmospheric structures could generate significant 
photospheric variability. Here, we estimate the level of thermal infrared vari- 
ability expected in successive secondary eclipse depths, according to hot Jupiter 
turbulent "shallow-layer" models. The variability, at the few percent level or 
more in models with strong enough winds, is within the reach of Spitzer mea- 
surements. Eclipse depth variability is thus a valuable tool to constrain the 
circulation regime and global wind speeds in hot Jupiter atmospheres. 



1. Introduction 



The regime of circulation in hot Jupiter atmospheres may be unlike any of the familiar 



cases in the Solar System. Hot Jupiters are gaseous giant planets found in c 
orbit s around their parent stars, with periods on the order of a few days (e.g. 



2006) . General arguments suggest that these planets should be tidally locked (jGuillot et al. 



ose, circular 



Butler et al. 



19961 ; iRasio et al.lll996l ; iLubow et al.lll997l ; lOgilvie fc Linll2004l ): their permanent day-sides 
are then continuously subject to intense stellar irradiation, while night-sides are only subject 
to modest internal energy fluxes. Given such an uneven energetic forcing, atmospheric 
winds would tend to redistribute heat around the planet. The nature and efficiency of this 
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redistribution p rocess is importan t in determining a variety of hot Jupite r observational 



properties (e.g. 


Seager et al. 


2005; 


Burrows et al. 


2005; 


Fortnev et al. 


2005: 


Barman et al. 


2005; 


Iro et al. 


2005; 


Burrows et al. 


12006: 


Fortnev et al. 


2006). 



Two groups have attempted to address the glob al atmospheric circulation problem in 



hot Jupiter atmospheres , using different approache s (IShowman fc Guillotl 120021 ; ICho et al. 



2003 



111 



Menou et al. 



200i ICooper fc Showmanl I2OO5I ) . On dynamical grounds, it was argued 



Menou et al.l (120031 . see al so Showman fc Guillot 2002), and explicitly shown via turbulent 



shallow-layer simulations in lCho et al.l (120031 . 120061 ). that tidally locked hot Jupiters occupy 
a regime of circulation that is qualitatively different from that of Solar System giant plan- 
ets. The few bands and prominent circumpolar vortices emerging in these simulations can 
be understood in terms of the large Rhi nes scale and Rossby de formation radius in these 
atmospheres (e.g., Cho fc Polvani 1996al lbl). iMenou et al. ( 2003 ) suggested that the large 
resulting spatial scales of moving circulation structures could lead to detectable hot Jupiter 
variability. 

Soon, interesting observational constraints will be placed on these circulation regime 
arguments. Three hot Jupiters have been detected through inf rared secondar y eclipses with 
the Spitzer Space Telescope: HD189733b ([Perning et al.l l2006l). HD209458b (foeming et al. 



20051 ). and TrES-1 (jCharbonneau et al.l 120051 ) . The planetary day-side flux is deduced from 
the eclipse depth measured when the planet is hidden behind its star. Repeated eclipse 
measurements could thus reveal detectable levels of variability of the planetary day-side flux 
in these three systems. In this Letter, we quantify the level of thermal infrared variability 



expec t ed in secondary eclipse depths, according to the shallow-layer models of ICho et al. 
(I2003l . l2006h . 



2. Shallow-Layer Circulation Models 



As explained in detail in lCho et al.l (120031 . 120061 ). turbulent equivalent-barotropic models 
published to date greatly emphasize dynamical aspects of hot Jupiter atmosph eric circula- 
tion. We adopt the same notation and adiabatic models as in ICho et al.l (120061 ): the global 
wind strength, U, and the net amplitude of radiative forcing, 77, are parameterized, not 
predicted. However, given these two (and other relevant global planetary) parameters, the 
turbulent atmospheric circulation is consistently found to develop a broad equatorial wind 
and two large circumpolar vortices revolving around the poles in several planetary days (= 
orbits). The dynamically modified layer thickness in the se models i s a pro xy for the planet's 
photospheric temperature field. We refer the reader to I Cho et al.l (120061 ) for details on the 
models, as well as a vast exploration of their parameter space. 
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Our limited goal here is to show that the thermal variability expected in at least some 
of these models is sufficiently large to be detectable via repeated Spitzer secondary eclipse 
measurements. Consequently, we focus on a limited set of four models, with global wind 
speed values from U = 100 to 800 m s _1 and a moderate amplitude of radiative forcing, 
7] = 0.05 (allowing the weak thermal contrast of features in slow wind models to remain 
apparent). 

The models are explicitly calculated for HD209458b parameters, at moderate T63 
(192 x 96 grid) resolution, over a hundred planetary days or more. Resolution tests (up to 
T170) show that T63 is sufficient to capture atmospheric temperature features well enough 
for our present purpose. Daily outputs from the simulations are used to generate model tem- 
perature maps of the day-side thermal emission in our variability study (see below). Table 1 
summarizes the range of photospheric temperatures derived from these four models, for the 
two brightest planets in our study. 

Figures [1] and [2] show snapshots of orthographically projected, day-side temperature 
maps in HD209458b models with U = 100 and 400 m s -1 , respectively, for two successive 
eclipses (i.e. after one HD209458b day). These projections illustrate how thermal variability 
in total eclipse depth is expected, from one eclipse to the next, if large, high contrast tem- 
perature features, associated with moving circulation structures, are present (in particular, 
the cyclonic circumpolar vortices most obviously visible in Fig. [2]). Each temperature map in 
Figs. [Hand [2] is shown partially eclipsed, for the specific geometry of the HD209458 system. 
According to these circulation models, cold polar vortices have relatively small (apparent) 
areas, so that the magnitude of their contribution to a variable eclipsed day-side flux is 
unclear without a detailed calculation. 



3. Thermal Variability in Eclipse Depth 



Our method to calculate ecl ipse depths in models like the ones shown in Figs. [T] and [2] 
follows very closely that used in iRauscher et al.l (120071 ) to study partial eclipse diagnostics. 
Planet-daily outputs from the circulation simulations provide the temperature fields used 
to model successive eclipses. Accounting for system specific inclination and geometry, these 
temperature fields are orthographically projected onto a 2D disk discretized with (100, 200) 
resolution elements in (r,9). To calculate spectra, we assume that the yertical temperatur e 
profile follows radiative equilibrium according to the cloudless models of ISeager et al.l (120051 ). 
under the assumption that the local flow temperature from the circulation model equals the 
effective temperature in the radiative model. We then integrate the spectral emission con- 
tributed by each apparent surface element on the planetary disk, in the global range of 
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effective temperatures from 700 to 2000 Monochromatic fluxes at Earth are then inte- 
grated over Spitzer spectral bands (for IRAC, IRS, and MIPS)0to predict the corresponding 
successive eclipse depths. In the future, it will be important to improve upon this simple 
treatment by fully incorporating radiative transfer in circulation models. 



As in iRauscher et al.l (120071 ). we also perform an idealized, bolometric blackbody anal- 
ysis, to avoid relying exclusively on model specific features of the cloudless spectral emission 
models used. In these simpler blackbody models, the bolometric flux contributed by each 
planetary disk surface element is scaled as the fourth power of the local photospheric tem- 
perature. 

Figure [3^ shows relative variations in successive eclipse depths predicted by the four 
HD209458b circulation models, assuming simple bolometric blackbody emission. Variations 
are semi-periodic, as expected from the quasi-periodic motion of the circumpolar vortices 
around the planet. As the magnitude of the global wind speed, U, is increased from 100 to 
800 m s _1 , leading to higher contrast motion- induced temperature structures, the amplitude 
of thermal infrared variability also increases, from a few % to ~ 60-70%. Figure [3b shows 
eclipse depth variations for the HD209458b model with U = 400 m s _1 as before, but this 
time using the detailed spectral emission models to calculate contributions in various Spitzer 
bands. While the overall variability scale is similar, it is clear that, by using Spitzer bands, 
one preferentially filters emission from a selective range of temperatures on the planetary 
disk, which contribute to varying levels of variability. The effect is substantial and we find 
that variability is the strongest in shortest wavelength Spitzer bands, where thermal contrast 
from the cold circumpolar vortices is the highest. 

We calculated explicit circulation models only for HD209458b but we can use the general 



dynamical similarity of hot Jupiter atmospheres (IMenou et al.l 120031 ) to rescale simply our 
results for HD189733b and TrES-1. Assuming identical dynamics but allowing for different 
average atmospheric temperatures, we linearly rescale our model temperature maps propor- 
tionally to the zero albedo, fully red istributed equilibrium temperature, T p , of the other two 



planets (as in IRauscher et al.l 120071 ). We find that the eclipse depth variability properties 
for HD189733b and TrES-1 do not deviate from those of HD209458b by more than ~ 10%. 
Increasing planetary albedos also has little effects on the variability properties, as long as 
albedos remain <C 1. Finally, as a matter of generality, we have checked that arbitrarily 
varying any system's orbital inclination in the range 80-90° has little effect on its variability 



1 For temperatures slightly below 700 K in the U =800 m s 1 models, a simple linear extrapolation of the 
spectra is performed. 

' http : //ssc . spitzer . caltech . edu/obs/ 
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properties. 



4. Detecting Eclipse Depth Variability 



We now address the feasibility of detecting variability in eclipse depth, at the level 
predicted by the above models, with Spitzer. Let us define a e d as the fractional error 
on the eclipse depth, that is the ratio of the full (1-cr) error on the eclipsed flux to the 
eclip sed flux itself. Values of a e d for the existing seco ndary eclipse measurements (taken 



from lCharbonneau et al.ll2005l : iDeming et al.ll2005l . 120061 ) are listed in bold in Table 2. These 
numbers already suggest that eclipse depth variability at the level of 5-20% could be detected 
in these systems. 

Our variability models indicate, however, that specific Spitzer bands may be much more 
useful than others for eclipse depth variability detections. We perform simple estimates of 
a e d errors for any combination of Spitzer instrument and hot Jupiter system as follows. We 
use the same set of system parameters as in Table 1 of iRauscher et al.l (120071 ). Assuming that 
errors on the non-eclipsed flux are comparatively very small, we write a e d = ai/VN, where 
o\ is the noise per data point in units of the eclipsed flux and N is the number of single data 
points collected during a full eclipse period. N is the ratio of the secondary eclipse duration 
to the instrumental cadence (taken from existing measurements). The eclipse duration time 
is calculated as t ec = 2a/ (R* + R p ) 2 — a 2 sin 2 (90° — i)/v, where v = 2na/P is the planet's 
orbital velocity, P its orbital period, a its orbital semi-major axis, i is the or bital inclination, 



and R * and R p are the stellar and planetary radii, respectively (see Table 1 of IRauscher et al. 
2007h . We obtain values of t ec = 1.76, 3.23 and 2.55 hrs for HD189733b, HD209458b and 
TrES-1, respectively. 

Finally, we extrapolate instrument-specific a e d errors known for one system to the other 
two systems of interest by assuming blackbody emission for both the star and the planet. 
This results in the following instrument-specific scaling between systems A and B: 



cr. 



ed 



a 



cd 




where d is the distance to the system and B\ is the Planck function evaluated at the central 
wavelength of the instrumental band under consideration. T* is the stellar effective tem- 
perature and T p = T^^J R*/2a is the fully redistributed planetary equilibrium temperature 



calculated in the small albedo limit, exactly as in iRauscher et al.l (120071 ). The resulting 
extrapolated values of a e d for the three systems of interest are listed in Table 2. 



Repeated eclipse depth measurements with IRAC for HD 189733b and HD209458b ap- 
pear to be the most likely to succeed in detecting atmospheric variability, at the few percent 
level or more. Such variability measurements would be difficult for TrES-1. Figure [3] shows 
that any pair of successive eclipses will generally not display the full range of variability 
amplitude: more than two eclipse measurements are needed to sample variability properties 
adequately. This requirement can be quantified by calculating distributions of fractional 
variations in eclipse depth over series of 2, 3, 4 or more successive eclipse measurements. 
Comparing these distributions of eclipse depth variations to the (l-a e d) errors listed in Ta- 
ble 2, we find that detecting the eclipse depth variability predicted by the U =400 or 800 
m s _1 models at the 2-3 a level requires 3-4 IRAC eclipses (at 4.5 or, slightly better, at 
8//m). A generally larger number of eclipses is needed to detect the variability predicted 
by models with smaller global wind speeds. Eclipse depth variations at the level predicted 
by the U =100 m s _1 model would be systematically masked by eclipse depth measurement 
uncertainties, according to our estimates. Finally, we note that a sufficiently large number 
of successive eclipse measurements could reveal the quasi-periodicity apparent in Fig. [3J 



5. Discussion and Conclusion 

We have illustrated, using simple turbulent shallow-layer circulation models, how eclipse 
depth variability can be used to constrain the circulation regime and global wind speeds in 
hot Jupiter atmospheres. 

Clearly, our circulation models and variability predictions are highly idealized. We 
have focused on simple thermal diagnostics in models describing an atmosphere as a single 



Iro et al. 


2005; 


Seaeer et al. 


2005; 



20061 ). the presence of high-altitude haze or the possible existence of clouds could all seriously 
affect our conclusions. Even strictly within the framework of our shallow-layer models, we 
know that the parameterized amplitude of net radiative forcing, r], can affect t he thermal 



contr ast, and therefore the detectability, of moving atmospheric structures (see ICho et al. 



20061 . for details). We have recalculated all our circulation and variability models with 
increased values of r) — 0.1 and 0.2jf| We find that the variability amplitude is reduced by 



3 In these models with stronger radiative forcing, an additional 5-10% contribution to the flow kinetic 
energy results from conversion of available potential energy. 
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up to a factor of two from the case with rj = 0.05 shown in Fig. [3j Values of rj > 0.2 would 
furth er reduce the variability amplitude. In the future, multi-wavelength phase curve data 
(e.g., lHarrington et al.l 120061 ) may provide useful observational diagnostics on the relevant 



value of rj for a given planet. 

Despite these shortcomings, our results are promising in showing that eclipse depth vari- 
ability is a new and potentially powerful tool to diagnose circulation and wind speeds in hot 
Jupiter atmospheres. In the future, as more refined atmospheric models are developed and 
more data becomes available, this tool should become increasingly useful in characterizing 
hot Jupiter atmospheres. 

We thank an anonymous referee for comments that helped improve the manuscript. 
This work was supported by NASA contract NNG06GF55G, NASA Astrobiology Institute 
contract NNA04CC09A and a Spitzer Theory grant, q 
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Table 1. 



Cho et al 



(120061 ) models under consideration 



Min-Max Temperature (K) 





U (m s- 1 ) 


HD209458b 


HD189733b 


Model 1 


100 


1147-1275 


956-1062 


Model 2 


200 


1112-1287 


926-1072 


Model 3 


400 


1006-1308 


838-1088 


Model 4 


800 


669-1372 


557-1141 
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Table 2. Measured (in bold) and estimated values of a ed . 



Instrument, wavelength 


HD 189733b 


HD 209458b 


TrES-1 


IRAC, 4.5 urn 


0.018 


0.039 


0.20 


IRAC, 8 fjm 


0.014 


0.026 


0.16 


IRS, 16 jan 


0.055 


0.11 


0.62 


MIPS, 24 jum 


0.086 


0.18 


0.96 
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U=100, first eclipse 



0=100, second eclipse 




Fig. 1. — Partially eclipsed temperature maps (in K) in a HD209458b model with a 100 
m s _1 global wind speed, for two successive eclipses. Little thermal infrared variability is 
expected. 




Fig. 2. — Same as Fig. [T] in a model with a 400 m s 1 global wind speed. Significant thermal 
infrared variability is expected in this case. 



13 




successive eclipses 



Fig. 3. — Relative flux variations at full eclipse according to HD209458b circulation models. 
Discrete points are connected for clarity and fluxes are normalized to the minimum value in 
each case, (a) Variations for simple bolometric blackbody emission, in models with [7=100 
(purple line), 200 (red), 400 (blue), and 800 m s -1 (black), (b) Variations according to 
detailed spectral emission models, in various Spitzer bands, for the [7=400 m s -1 model of 
HD209458b. From top to bottom, the curves correspond to 3.6/im (black), 4.5/mi (blue), 
8/im (red), 6/mi (green), 16/im (orange), and 24/mi (yellow). 



